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Abstract 

The structure of Pais-Uhlenbeck oscillator in the equal-frequency limit has 
been recently studied by Mannheim and Davidson [Phys.Rev. A71 (2005), 042110] 
It appears that taking this limit , as presented in the above paper, is quite subtle 
and the resulting structure of space of states - involved. In order to clarify the 
situation we present here the proper way of taking the equal- frequency limit, 
first under the assumption that the scalar product in the space of states is pos- 
itive defined. We discuss also the case of indefinite metric space of states. We 
show that, irrespective of the way the limit is defined, the limiting theory can be 
hardly viewed as satisfactory. 
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Pais-Uhlenbeck quartic oscillator [T] is described by the Lagrangian 

m , 2 mu 2 2 mX.. 2 
L = -q-—q-—q (1) 

Its behaviour depends on actual values of parameters to and A. In what follows we keep 
m and u fixed while varying A. The relevant equation of motion reads 

with 



2 1 ± VI - 4A^ 2 
"i, 2 2A (3) 

For large A (A > ^-y) both frequencies are complex. On the other hand, in the range 
< A < they are real; in the limiting case A = there is a double degeneracy 
uj\ = uj 2 = V2co. Finally, if A < 0, one frequency is real while the second one - purely 
imaginary. 

In order to quantize our theory one has to put it first in Hamiltonian form. This 
can be achieved within Ostrogradski formalism [2], [3]. It is well known [5], [6] 
that the Ostrogradski procedure is essentially a form of Dirac method for constrained 
theories. 

In our case one finds [Tj, [8] -j- [ID] the following canonical variables 

qi = q, q2 = q 
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= —mXq 




dq 



Ili = — = — - — ( — ) = m(q + A q) 

(4) 



together with the Hamiltonian 



,2 , mUj2 2 '" 2 



Quantization can be now performed in the standard way by imposing the commutation 
rule 

[qi,Hj] = ihSij (6) 

Consider first the range < A < then u\ > uj\ > 0. To make the structure 
of the Hamiltonian transparent we perform the following canonical transformation (cf. 
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Ref.pQ) 



qi = , =(-xi + x 2 ) 



n 2 = m\\— 2 2 (— + W2X2] 



Note that the above transformation becomes singular in the doubly degenerate limit 
A — > jij. In term of new variables the Hamiltonian takes particularly simple form 

The eigenvectors of H are uniquely determined (up to a phase factor) by two nonneg- 
ative integers ri\ , n 2 

H I ni,n 2 ) = (-%ui(ni + i) + ^a; 2 (n 2 + ^ | ni,n 2 ) (9) 

The spectrum of i7 is simple provided ^ is irrational; for rational ^ (superintegrable 
case) there is a degeneracy. 

The wave functions in the coordinate representation read 

(x 1 ,x 2 I ni,n 2 ) = (10) 

where N(n) = (y/n ■ 2 n ■ n\)^^ . 

The spectrum of H, as given by eq.Q, is unbounded from below. One gets positive 
energy spectrum by admitting indefinite metric in the space of states. To this end we 
consider the space of states endowed with the positive-definite scalar product (-, •) and 
define the "physical" scalar product with the help of metric operator 77 

($ I = ($, 77*), r] = r] + = r]' 1 (11) 

Denoting by " *" the hermitean conjugate with respect to the scalar product (■ | ■) one 
finds for any operator A. 

A* = r]A + rj (12) 
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Let ai,af be the creation/anihilation operators constructed out of x^pi 



Xi 1\ 



h , 

- \ CLi ~ 



2mw 



We define 



Then 



P^J^i^af) (13) 



r) = (-1)^ = e ^ ai (14) 



(ni,n 2 | n'iX 2 ) = {-^) ni ^nxn'^n 2 n' 2 \ (15) 

moreover, x* = Xi, p* = pi imply 

xf = {-!)%, pt = {-±)% (16) 

Defining (c.f. pQ) 

X — I 7, .'77 ~^ j ^ — 3/ 2 

Pi = =F^i, P 2 =P2 (17) 
we find that x'^p^ are hermitean (with respect to "+" conjugation) and 

M£+^ f MS +! ^ ? ) (18) 

Therefore, the spectrum of H is now positive definite 

E num = hujtini + -) +huj 2 (n 2 + ^) (19) 

The "physical" subspace is spanned by the vectors | 2ni,n 2 ). 

Let us now consider the degenerate case A = To reveal the structure of the 
Hamiltonian we define new variables Qi , A by [1] 

Qi A 

Qi = — ?= H 

2\/2 rau; 

92 = -r- + (20) 

2 m 

fr ^ 3 n 
III = - -mojQ 2 

lt 2 = i=mQi 

2u 4^2 ^ 



Eq. (120]) defines a A independent canonical transformation. In terms of new variables 
the Hamiltonian takes a particularly simple form 
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H = y/2u{Q x P 2 - Q 2 A) - + O2) ( 21 ) 

Again, it is a sum of two commutating pieces: the first proportional to the angular 

momentum operator while the second represents the length of Q squared. Therefore, 
the energy is the sum of discrete and continuous parts and depends on the combination 
of two quantum numbers [1] . 

In what follows we find it convenient to use the momentum representation: Pj — > 
Pi, Qi ~ ihijp:- Let us define the polar coordinates in momentum space by P\ = 
PcosO, P2 = PsinB. Then 

■ r , 9 muj 2 h 2 ( d 2 Id 1 d 2 \ 
H = + — + + j ; (22) 

the spectrum and normalized eigenf unctions read 

„ . ( 1- mu>hk 2 \ , . 

E nM = uhW2n ^— j (23) 

*n,k( p ) = ^ j ^ kp y n& 

We shall now consider the limit A — ► p. This is slightly subtle due to the fact that 
the spectrum of the Hamiltonian changes in this limit from discrete into continuous 
one. Let us put 

1 - 4Xu 2 = e 2 , 0+ (24) 

Then 

V2u (l ± (25) 



^1,2 



The energy spectrum (j5J) can be rewritten as 



E nun2 = V2uh(n 2 - nx) —u)he(ni + n 2 + 1) (26) 



In the limit e — > + the energy seems to be given by the formula [TT] 

E nun2 = V2uh(n 2 - fix) (27) 

This is, however, not the case. We should take into account that the limiting Hamilto- 
nian has a continuous spectrum given by eq.( f24|) . Therefore, the proper way of taking 
the limit is to let n\ , n 2 — > 00 in such a way that 

n = n 2 — Hi 

muhk 2 , . 

e{n 1 + n 2 ) = —j=- (28) 



are fixed. 

We shall show that, indeed, by considering this limiting procedure one recovers the 
wave functions ( 1231) of degenerate Hamiltonian. 

Let us note that we cannot take the equal frequency limit directly. This is due to 
the fact that the very coordinate representation becomes singular in this limit as is 
clearly seen from eq.fH). On the other hand, the momentum representation based on 
Pi,P 2 (cf- eqs.(l22l. (|23|) ) is always well-defined. Therefore, the first step will be to 
rewrite our wavefunctions ( JTOl) in momentum representation. To this end we write Pi 
and P2 in terms of x^s and p' { s. By virtue of eqs.([7]) and (1201 we find 



Pi 



P> 



b 2 pi + bip 2 



\/2co(—biXi + 62^2) \/2u 

2V2 + A "* 2 ) 



-ih ( d , d 
b 2 ^ \- h 



dxi dx 2 J 
-61X1 + b 2 x 2 ) 



(29) 



m /3 
61 2 = — = I 



It is now straightforward to find the relevant transition functions (x\,x 2 \ P\,P 2 ) by 
solving the corresponding eigenvalue equations 



Pi(xx,x 2 I Pi,P 2 ) = Pi{x x ,x 2 I Pi,P 2 ) 
Using eqs.( l29i) one obtains 
(x 1 ,x 2 I Pi,P 2 ) = 



(30) 



(31) 



m\(ujj — ul) 



5 -b 1 x 1 + b 2 x 2 



A(w? - u%)P 2 



" \ 2V2nhu " I ^ ' x/2oo 
This allows us to pass to the momentum representation 

(Pi,P 2 I nx,n 2 ) = J dxidx 2 (P u P 2 | x 1 ,x 2 )(x 1 ,x 2 | ni,n 2 ) 



(32) 



Doing one integration with the help of delta function one arrives at the following rather 
complicated expression: 



(Pi, P 2 I n u n 2 ) 



N{n\)N{n 2 ) Jm 2 ^^ 



bl + bl 



m\(u 2 



P 2 +iPl 



h 2 \ 2x/2nhu 



(i,2 + f,|)2 2 



dvH ni 



muji 
~1T 



b 2 v - bi 



V 



b 2 + b 2 



muj2 



(33) 
\ 



b 2 + b 2 



X(u 1 u> 2 ) I b 1 b 2 (ui 2 -^i) 



P2+iPl 
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In principle, the last integral could be also taken. However, this is not necessary. We 
can take the limit e — > + directly in eq. (1331 . In this limit 

»u - ^ (l ± j) (34) 



(35) 



Keeping dominant terms in eq. (l33l) one finds 

N( ni )N(n 2 )^ 




: 'P 2 + P x ))H n2 [y-J - J - (-P 2 + iPi) e- 



2y/2mhu J 2 \ y 2\]2mhio 

The last integral is taken explicitly yielding [12 



(Pi,P 2 | na,n 2 ) ~ (36) 

n 



(for definiteness we have assumed here n\ > n 2 ; the opposite case goes along the same 
way). Now, it is easy to take the limit e — > 0, — > oo with n = n 2 — n\ and 
e{n\ + n 2 ) = fixed. Using Stirling formula and the asymptotic form of Laguerre 

polynomials [12J we arrive finally at the following result 




By comparying eq. (!36l) and (|23~1) we conclude that by taking the limit of equal frequen- 
cies in the way prescribed above we recover the wavefunctions of degenerate Hamilto- 
nian. To complete the arguments let us only note that the additional y/e factor comes 
from the fact that the discrete eigenfunctions have unit norm while the norm of limiting 
one is infinite. 

Let us now consider the question whether the double-frequency case can be quan- 
tized in the way which yields positive-definite Hamiltonian at the expense of having 
indefinite metric in the space of states. More precisely, the problem is the following: for 
differening frequencies we choose the positive-energy quantization scheme; is it possible 
to perform the equal-frequency limit keeping the energy positive? 

In order to answer this question one has to adopt some definition of taking the 
equal-frequency limit. The most reasonable way is to keep the initial operators ^ , 11^ 
as given once forever linear operators acting in some linear space of states and varying 
the Hamiltonian by varying some of its parameters (A in our case). For e > one can 
define, via eqs.flT]) and ([TBI) , the hermitean conjugation " +" for q^s and U^s (and, hence, 
convert the space of states into usual Hilbert space) in such a way that the Hamiltonian 
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(jSJ) becomes a positive-definite operator. The definition of hermitean conjugation " +" 
depends, of course, on e. It is easy to compute, using ([7]) and f|T6|) . that, for example 

qt ^-Jqi- (38) 



e 2 v m 

Therefore, the reasonable conjugation rule cannot be imposed in the limiting case. One 
can support this conclusion by considering the operators a and b defined by eqs. (33) 
of Ref. [8]. In fact, it is easy to check that a diverges if the limit is taken in the way 
prescribed above; this is not in contradiction with the fact that their commutators and 
the Hamiltonian (eqs. (35), (40) of Ref [8J) are well-behaving in the limit e — > 0. 

However, one can argue that our prescription of taking the e — > limit is not the 
only possible. In fact, one can argue that it is sufficient to get the regular limiting 
commutation rules and the Hamiltonian. This seems to be the strategy adopted by 
the authors of Ref. [8] . 

Let us consider in some detail their construction. The limiting commutation rules 
and the Hamiltonian read 

[a, a*} = 0, [b, b*} = 0, [b, a*] = /j, [a, b*} = ^ [a, b] = 

H — —(2b*b + a*b + b*a) + const (39) 



By defining 



A l = -j={a + b), A 2 = -==(a-b) (40) 



one obtains 

[A, Aj] = 0, [A$,A$=0, [A,A*} = 0,i^j, [^,^ = 1, [A 2t A$ = -l 
H = uj{2A\A x - A\A 2 - A\A X ) (41) 

Then we obtain the indefinite Fock space for two degrees of freedom. Define new metric 
operator 

t = e mA i M (42) 

and new conjugation 

B + = tB*t (43) 

Then Af = A*, A\ = —A\ and one obtains standard commutation rules together 
with Fock representation with positive metric and " + " playing the role of hermitean 
conjugation. 

Now, H can be viewed as acting in standard Hilbert space, 

H = u(2A+A 1 - A\A 2 + A+A-t) (44) 
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We see that H is not normal 

[H,H+}^0 (45) 

so it cannot be diagonalized. Moreover,, 

[H, N] = (46) 

where N = AfA 1 + A 2 A 2 is the total number operator. Therefore, we can consider 
H as acting separately in each finitedimensional subspace of fixed eigenvalue of N . It 
is represented as a matrix in any such subspace and we want to determine its Jordan 
form. To this end consider the subspace spanned by | ni, n 2 ), n\ + n 2 = n; then 

(H - un) n+l I m, n 2 ) = {H- uN) n+1 \ n u n 2 ) = 

= u n+ \A+A 1 - A+A 2 - A+A 2 + AtA x ) n+1 \ n x n 2 ) = 

= u n+1 ((At + At){A, - A 2 )) n+l I m, n 2 ) = 

= u n+l (Af + At) n+ \A 1 - A 2 ) n+1 I m,n 2 ) = (47) 

because [Af + A 2 , A\ — A 2 ] = 0. Therefore, in any subspace corresponding to the 
eigenvalue n of N H takes the form of single Jordan block; any such subspace contains 
exactly one eigenvector corresponding to the eigenvalue nuj. In order to construct it 
we note that 



lH,A+ + Aj]=u(A+ + A+) (48) 

Therefore 

\n) = (At + A+) n I 0,0) (49) 
is the eigenvector of H corresponding to the eigenvalue ntu of H. Now, using 

I "> = p Q VUy/(n-k)\ I k,n- k) (50) 

one finds 

(n\n)=n\±(^j(-lT~ k = (51) 

We see that all eigenvectors of H have zero norm. However, it is reasonable to assume 
that the physical space of states is spanned by the eigenvectors of H. This is a subspace 
of zero norm so it can hardly be viewed as physical one. 

Concluding, we see no way to get the positive-energy quantized degenerate PU 
theory, at least by limiting procedure. 
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